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ABSTRACT 


The  effective  interaction  in  nuclear  matter  is  defined  as  the 
long  range  part  of  the  two  body  potential.,  which  in  the  Born 
approximation  gives  the  single  particle  potential  energy  for  the 
average  momentum  in  the  Fermi  sea.  For  the  Gammel  -  Thaler 
potential  the  effective  interaction  has  been  calculated,  first 
for  the  free  particle  propagator  and  then  for  the  nuclear  spectrum 
The  result  shows  that  in  the  free  particle  case  the  separation 
distance  £0  is  constant  over  a  wide  range  of  densities  and  does 
not  lead  to  saturation.  The  nuclear  separation  distance  £(Kp) 
changes  quite  rapidly  with  the  Fermi  momentum  Kp,  and  while  for 
low  densities  it  is  very  close  to  for  higher  densities  it 

becomes  very  much  larger.  At  the  density  corresponding  to 
Kp  =  1.5F“1,  the  long  range  potential  starts  at  ^=1.166f-1. 

The  minimum  of  the  total  energy  per  particle  occurs  at  Kp  =  1.43F” 
and  is  about  -lOMeV. 

A  variational  method  is  also  used  to  calculate  the  reaction 
amplitude  for  nuclear  matter.  It  is  found  that  the  potential 
energy  for  state  is  -15.45  MeV  as  compared  to  -l4.  MeV 

obtained  by  perturbation. 
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Chapter  I. 

INTRODUCTION 

The  study  of  the  properties  of  nuclei  in  its  most 
general  terms  includes  a  large  variety  of  topics  ranging 
from  the  study  of  nuclei  in  their  ground  states  to  scattering 
and  reaction  processes  at  high  energies.  In  this  paper  we 
shall  concern  ourselves  with  a  much  more  restricted  aspect 
of  nuclei,  the  binding  energy  of  nuclear  matter. 

The  semi-empirical  mass  formula  developed  by 
Bethe  and  Weizs&cker  around  1935-36  describes  the  "averaged 
binding  energy  curve  of  atomic  nuclei  with  atomic  mass  A 
and  atomic  number  Z." 

p  /  p  P—1  P-1 

Binding  Energy  =  -a^A  +  a^A  +  a^Z  A  +  1/4  a^(A/2  -  Z)  A 

El  -  (  8MeV)  A 

where  a^  is  a  measure  of  the  true  nuclear  interaction 
energy,  a^  measures  the  "unsaturated"  nuclear  binding 
energy  at  the  nuclear  surface,  a^  is  determined  by  the 


_ 


2. 


Coulomb  energy  between  protons  and  a,,  is  a  repulsive 


term  which  says  that  the  most  stable  nuclei  have  Z  =  A/2. 


Green 


(1) 


has  made  the  most  accurate  fit  to  the  semi- 


empirical  mass  formula.  In  particular,  he  has  found 
the  value  a^  =  15. 5  MeV. 

The  concept  of  nuclear  matter  is  an  idealization 
based  on  the  notion  that  a  large  number  of  nucleons 
would  form  a  stable  configuration  if  Coulomb  repulsion 
between  protons  did  not  exist.  The  properties  of  this 
configuration  will  be  independent  of  the  number  of  nucleons 
if  that  number  is  so  large  that  surface  effects  can  be 
neglected.  From  the  semi-empirical  formula  it  is  easy 
to  see  that  the  binding  energy  per  nucleon  is  equal 
to  -a^  because  for  stable  nuclei  with  large  number  of 
nucleons  all  of  the  terms  except  the  first  term  can  be 
neglected . 


Nuclear  matter  may  be  looked  upon  as  a  gas  of 
strongly  interacting  Fermi  particles,  and  an  individual 
particle  or  shell  model  of  heavy  nuclei  consists  of  a 
description  in  terms  of  momentum  states  with  an  altered 
energy  -  momentum  relation.  Ordinary  perturbation  theory, 
with  free  particles  as  the  unperturbed  system,  serves  as 
a  crude  starting  point  by  which  one  might  attempt  to 
establish  a  theoretical  basis  for  the  shell  model. 

In  such  a  treatment,  we  would  take  the  entire  interaction 


' 


- 
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as  a  perturbation.  This  was  done  by  Euler  as  early  as 
1937,  and  later  by  Huby .  The  results  were  discouraging 
because  the  calculated  binding  energy  was  too  small  and 
the  convergence  rate  was  poor. 

Three  more  successful  independent  calculations 
of  the  binding  energy  of  nuclear  matter  have  been  made 

( p ) 

with  realistic  nuclear  forces,  by  Brueckner  and  Masterson 

( 8 ) 

by  Razavy  following  the  reference  spectrum  method 

( 4 ) 

developed  by  Bethe,  Brandow  and  Petschek  and  by  Scott 

(5) 

and  Moszkowski 

The  potential  energy  of  the  individual  particles, 

at  least  sufficient  for  finding  the  total  nuclear  energy, 

might  be  determined  by  a  Hartree-Pock  self-consistent 

field  method.  However,  the  evident  inadequacy  of 

ordinary  perturbation  theory  and  the  nature  of  the  nuclear 

forces  (hard  cores,  exchange  character,  etc.)  leads  one 

to  believe  that  in  any  accurate  calculation  one  must 

include  some  of  the  strong  correlations  between  the 

( 8 ) 

particles.  Brueckner  and  co-workers  have  developed 
a  method  of  treating  the  many  body  problem  which  takes 
a  first  step  in  this  direction. 

The  theory  of  Brueckner  and  co-workers  permits, 
in  principle,  calculation  of  properties  of  complex  nuclei 
in  terms  of  the  potential  between  two  nucleons.  The 
theoretical  foundation  of  the  theory  has  been  given  by 


' 

as 
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Goldstone ^ ^ and  by  Hugenholt z ^  ^ .  Goldstone  uses  an 
expansion  of  the  development  operator  whereas  Hugenholtz 
expands  the  resolvent  operator.  Brueckner  and  Gammel^^, 
aside  from  further  developing  the  method  did  extensive 
numerical  calculations,  using  the  best  results  then 
available  in  the  interaction  between  two  nucleons,  and 

obtained  results  in  good  agreement  with  experiment. 

( 2 ) 

Brueckner  and  Masterson  used  an  approximation  of 
the  Brueckner  formalism  which  is  simpler  than  that  of 
Brueckner  -  Goldstone. 

The  Brueckner-Goldst one  (BG)  method  is 
motivated  by  the  fact  that  the  nuclear  forces  turn 
strongly  repulsive  at  short  distances.  This  makes  any 
kind  of  perturbation  treatment  inapplicable,  since  the 
matrix  elements  of  an  infinite  repulsive  core  potential 
with  respect  to  independent  particle  wave  functions  are 
all  infinite.  However,  it  is  well  known  that  a  potential 
containing  a  hard  core  does  not  have  an  infinite  effect 
on  the  wave  function  of  the  interacting  pair.  In 
particular,  we  obtain  a  finite  phase  shift .  The  two- 
particle  wave  function  is  simply  "pushed  away"  from  the 
region  of  the  core. 

Brueckner' s  theory  assumes  that  while  the 
correlations  between  any  two  particles  can  be  quite  large, 
any  correlations  involving  more  than  two  particles  at 


5. 


a  time  can  be  neglected.  This  seems  to  be  a  quite 
good  approximation  for  nuclear  matter.  Accordingly, 
the  nuclear  potential  energy  is  just  the  sum  of  the 
interaction  energies  of  each  pair  of  nucleons  and  the 
latter  is  determined  in  terms  of  the  two-particle 
scattering  properties. 

The  Brueckner’s  theory  can  be  criticized 
from  two  points  of  view:  From  the  basic  point  of  view, 
it  may  be  questioned  whether  it  actually  leads  to  ground 
state  of  nuclear  matter,  and  from  the  practical  point 
of  view  the  numerical  calculations  required  are  complicated 
and  not  very  transparent.  On  the  problem  of  the  basic 
side,  Luttinger,  Kbhn  and  Ward ^ 10 5 11 ^  have  shown  that 
the  BG  method  is  indeed  correct  to  all  orders  of 
perturbation  theory  if  (1)  the  particles  have  spin  1/2, 

(2)  the  interaction  between  them  is  isotropic,  and  (3) 
the  Fermi  surface  without  interaction  is  isotropic  in 
momentum  space.  These  conditions  are  evidently  satisfied 
for  nucleons  interacting  with  central  forces.  For  the 
case  of  nucleons  interacting  with  tensor  forces,  they 
showed  that  the  BG  theory  is  correct  in  the  second  order. 
Another  problem  on  the  basic  side  is  whether  perturbation 
theory  is  valid.  The  chief  argument  of  it  is  that  nuclear 
matter  may  exhibit  a  phenomenon  similar  to  super- 

(  4 ) 

conductivity.  Bethev  has  pointed  out  that  the  "super¬ 
conductivity  pairing"  has  a  negligible  influence  on 


. 


the  average  binding  energy  of  nuclei,  so  that  the  BG 

method  can  be  considered  as  sufficiently  established 

on  theoretical  grounds.  However,  the  numerical  work 

is  complicated,  especially  due  to  the  repulsive  core  in 

the  nucleon  interaction.  Levinger,  Peierls,  Razavy  and 

(  ir) 

collaborators v  ^'have  therefore  suggested  to  ignore 
the  repulsive  core  and  to  replace  its  effect  by  a  velocity 
dependent  potential.  There  is  no  objection,  in  principle, 
against  the  use  of  velocity-dependent  potentials.  Actually 
a  priori  it  is  hard  to  tell  whether  a  velocity-dependent 
or  a  static  potential  gives  better  results.  But  Charap, 


Fubini,  and  Tausner 


(14) 


have  shown  that  at  least  at 


moderate  energies,  dispersion  theory  leads  to  an  un¬ 
ambiguous  definition  of  a  static  potential  which  can  be 
calculated  in  terms  of  exchange  of  one,  two,  etc.,  pions 
(see  however  Ref.  15).  In  nuclear  matter,  the  filled 
states  have  momenta  up  to  =  1.5  which  is  moderate. 

So  that  an  ordinary,  static,  nucleon-nucleon  potential 
should  give  a  good  description  of  nuclear  matter. 

A  great  simplification  over  the  direct  solution 
of  Brueckner  integral  equation  was  made  by  Moszkowski 
and  Scott  ^  (MS).  They  developed  an  approximation 
scheme  which  makes  use  of  some  of  the  physical  features 
of  the  nuclear  forces.  The  essential  thing  is  that  the 
nuclear  forces  are  basically  attractive  at  moderately  low 
energies.  Although  they  apparently  turn  strongly  re¬ 
pulsive  at  short  distances,  the  attraction  at  larger 
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distances  more  than  compensates  this  repulsion.  In  other 
words,  the  repulsion  effectively  cancels  a  part,  but  not 
all,  of  the  attraction. 

The  nucleon-nucleon  interaction  is  therefore 

separated  into  a  short-range  part,  Vs  and  a  long-range 

£, 

part  V  .  The  short-range  part  is  defined  such  that  it 
gives  zero  phase  shift  for  free  particle  scattering  and 
is  represented  by  a  reaction  matrix  Gp  which  is  not 
very  different  for  nuclear  matter  and  for  free  nucleons. 

The  long-range  part  V  is  the  remaining  part  of  the 
interaction  which  can  be  treated  by  Born  approximation 
in  nuclear  matter  while  for  free  nucleons  due  to  the  Pauli 
effect  this  is  not  possible.  The  separation  between  short 
and  long  range  is  made  in  such  a  way  that  the  short  range 

g 

reaction  matrix  for  free  nucleons,  Gp, vanishes.  Then  the 
total  reaction  matrix  G  consists  mainly  the  first  and  second 
Born  approximation  for  V^,  the  difference  between  G?  in 
nuclear  matter,  G^,  and  for  free  nucleons,  Gp  ;  in  addition 
there  are  some  small  corrections. 

This  method  is  much  simpler  than  the  direct 
solution  of  the  reaction  matrix  of  Brueckner-Goldstone , 
but  Kfthler ^ ^ ^  finds  that  the  MS  calculation  of  G^  -  Gp 
may  not  be  sufficiently  accurate.  (Gjjj  is  the  nuclear 
reaction  matrix.) 

The  MS  method  has  been  improved  by  Bethe  et  al 
(B.B.P.)^).  They  also  have  developed  a  method,  the 


.  ■  - 
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reference  spectrum  method,  which  is  even  simpler  than 
the  MS  method.  They  replace  the  single  particle  energy 
for  the  intermediate  states  by  reference  energy  defined 
as  E^(k’)  =  A  +  k'^/2m*  where  m*  is  an  effective  mass 
and  k’  is  the  single  particle  momentum.  With  this 
replacement  the  energy  denominator  will  be  a  simple 
quadratic  function  of  the  relative  momentum  of  the  ex¬ 
cited  states.  In  addition,  if  the  Pauli  operator  is  also 
dropped  the  integral  equation  can  be  reduced  to  a  second 
order  differential  equation  which  can  be  handled  quite 
easily.  The  reference  spectrum  parameters,  A  and  m* ,  are 
chosen  to  fit  the  actual  particle  energies  with  the 
reference  spectrum  as  accurate  as  possible  for  k1  >  k^. 

The  reference  reaction  matrix  G  is  then  used  to  calculate 
the  Brueckner  matrix  G.  The  difference  G  -  G  is  due 
mainly  to  the  effect  of  the  exclusion.  The  spectral 

correction,  which  arises  from  the  approximation  of  the 
energy  denominator  is  quite  small  and  shows  that  a  detailed 
knowledge  of  the  form  of  the  single  particle  energy  is  not 
important  for  the  calculation  of  the  reaction  matrix.  The 
actual  energy  spectrum  must  of  course  be  calculated  from 
the  reaction  matrix. 

The  MS  separation  is  useful  to  assess  the  sensi¬ 
tivity  of  the  nuclear  binding  energy  to  the  behaviour  of  the 
nucleon  interaction  at  short  distances.  In  the  MS  method, 


.  't 
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the  short  range  forces  are  mainly  important  for  the  deter¬ 
mination  of  the  separation  distance  d.  We  can  also  deter¬ 
mine  d,  for  a  given  relative  momentum  k  ,  from  the  observed 
phase  shift  of  nucleon-nucleon  scattering  for  momentum  k  , 
if  we  assume  that  the  long-range  nuclear  potential  (say  for 
distances  greater  than  IF)  is  known.  We  can  simply  in¬ 
tegrate  the  Schrftdinger  equation  inwards  from  infinity 
(or  some  large  value  of  r),  until  we  reach  the  point  at 
which  the  logarithmic  derivative  of  the  wave  function  is 
equal  to  that  of  free  nucleons;  this  is  then  the  MS  separation 
distance  for  non-interacting  nucleons.,  dp. 

Puff  and  Martin  and  Mohling^^)  haVe  proposed 

a  new  method.  In  their  method  in  the  scattering  into  inter¬ 
mediate  states  the  exclusion  principle  is  not  taken  into 

account  and  the  excited  state  energies  are  replaced  by 
o 

k  /2M.  The  energies  of  holes  are  determined  self -consist¬ 
ently  as  in  other  theories.  The  Puff -Martin  theory  in¬ 
volves  partial  summations  and  cannot  be  developed  in  per¬ 
turbation  theory,,  even  when  well  behaved  potentials  are 
employed.  This  is  because  it  leads  from  zero-order  continuum 
states  to  a  resulting  state  which  is  bound.  It  is  well 
known  that  this  cannot  occur  in  perturbation  theory  taken 
to  any  finite  order.  On  the  other  hand,  the  Brueckner 
theory  can  be  expanded  in  perturbation  theory  when  weak 
well-behaved  potentials  are  employed. 


•  •i 
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The  main  task  of  this  paper  is  to  develop  an  easier 
method  to  calculate  the  Brueckner  reaction  matrix  which 
minimizes  the  effect  of  the  purely  phenomenological  part 
of  the  potential.  As  in  the  MS  method,  we  separate  the 
two  body  interaction  into  a  short-range  and  a  long-range 
part.  But  we  modify  the  definition  of  separation  point 
and  treat  the  Hamiltonian  of  the  short-range  interaction 
as  a  perturbation.  The  separation  distance  is  determined 
by  the  requirement  that  the  first  order  contribution  of  the 
short-range  potential  to  the  reaction  matrix  be  zero. 

This  definition  of  separation  distance  differs  from  that  of 
MS.  With  our  separation  the  first  order  expression  for 
the  reaction  matrix  is  just  the  first  Born  approximation 
for  the  long-range  part  of  the  potential  and  the  energy 
of  the  system  can  be  obtained  from  the  effective  Hamiltonian 
which  contains  the  long-range  interaction  and  the  kinetic 
energy.  The  details  of  solving  the  first  order  Brueckner 
equation  are  contained  in  Chapter  2.  The  first  order 
correction  is  mainly  due  to  the  reaction  matrix  of  the  short- 
range  part  treated  as  perturbation  which  is  discussed  in 
Chapter  3.  In  the  same  chapter  the  second  order  correction 
is  also  discussed.  Chapter  4  contains  the  numerical 
calculation  and  comparison  of  the  results  with  other  works. 


.  ,4-  v  -  f 
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We  also  use  a  variational  principle  for 

calculating  the  reaction  matrix  based  on  Schwinger’s 

method.  Hulthen^2^  has  proposed  a  workable  variational 

('21')  (22) 

method  which  has  been  modified  by  K&hn  .  Schwinger 
approached  the  problem  in  a  different  way.  While  both 
Hulthen  and  K8hn  used  the  differential  form  of  the 
Schrtidinger  equation,  Schwinger  started  out  by  taking 
the  appropriate  integral  equation  corresponding  to  the 
Schrbdinger  equation.  This  gave  rise  to  a  complicated 
functional  involving  Green’s  function. 

In  order  to  apply  the  Schwinger’s  method  in 
nuclear  matter  problem  we  use  the  integral  equation  for  the 
reaction  matrix  but  with  Greenes  function  different  from 
that  of  Schwinger’s  to  calculate  the  reaction  amplitude. 
This  variational  method  is  discussed  in  Chapter  5  and  the 
result  is  compared  with  that  of  perturbation. 


12. 


Chapter  2. 


THE  FIRST  ORDER  APPROXIMATION 


The  basic  assumption  in  this  theory  is  that  the 


correlations  between  any  two  particles  can  be  very  large, 
while  correlations  involving  more  than  two  particles  at  a 
time  are  small  and  to  the  first  approximation  can  be  neglected. 
Thus  the  total  Hamiltonian  of  the  nuclear  matter  contains 
the  sum  of  kinetic  energy  of  the  particles  and  two-body 
interactions  and  can  be  written  as: 


H  =  ?  T.  +  1/2) 

ii  ^ 


(2.1) 


In  equation  (2.1)  A  is  the  number  of  nucleons,  T^  the 
kinetic  energy  operator  for  the  particle  in  the  state  i 


and  V(r..)  the  potential  between  two  particles  in  state  i 


and  j  respectively.  The  interaction  can  be  divided  into 
two  parts,  i.e.  the  long-range  attractive  part  and  the 
short-range  part  which  is  strongly  repulsive  for  very  small 
r  and  then  strongly  attractive.  The  separation  distance 
which  is  not  yet  specified  is  denoted  by  E. . .  The 
Hamiltonian  can  be  written  as: 


H  =  H  +  H’ 

o 
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with 


H  =  I  T.  +  1/2  Z  V(r .  .)  e  (r,  . 
0  i  1  U  1J  10 


H'  1/2K  V(rlj)  9  (Clj  rij  ^ 


(2.2) 


Ho  is  unperturbed  Hamiltonian  and  H'  perturbation. 

The  total  Hamiltonian  H  defines  the  exact  wave 
vector  |  >  and  the  total  energy  E  of  the  ground  state 

of  the  system 


H  I  ip  >  =  E  I  ip  > 
1  o  1  r  o 


(2.3) 


and  Hq  defines  a  set  of  unperturbed  wave  vectors  |  > 


with  eigenvalue  E 


n 


H  $  >  =  E  $  >  . 

o  1  n  n  1  n 


(2.4) 


We  normalize  I  u>  >  such  that  <  $  I  i|<  >  =  1.  We  can 

1  0  0  0 

expand  |  ij;  >  in  terms  of  the  complete  orthonormal  set 


$  >  to  get 

n 


y  >  *  $  >  + 

o  1  o 


z: 

n  ^  o 


$  >  <  $  H '  ¥  > 

n  n  1  1  o 


E  -  E 


n 


(2.6) 


and  the  total  energy 


E  =  E  +  <$  H’  Y  > 
o  o  1  1  o 


=  E  +  AE. 

o 


(2.7) 


The  energy  shift  AE  is  a  function  of  the  separation 


distance  £>±.  .  If  =  c  where  c  is  the  radius  of  the  hard 
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core  (region  of  strong  repulsion),  then  AE  >  0 .  On  the 

other  hand  if  is  very  large  then  for  the  existence  of 

the  bound  state  AE  should  be  negative.  Thus  it  is  possible 

to  find  a  point  E, .  .  in  such  a  way  that  the  energy  shift 

J 

vanishes . 


The  wave  vector  I  i . . >  for  each  pair  of  particles 

1  ij 


is  defined  such  that 
pairs  except  ij 


V . . >  contains  correlations  of  all 

ij 


ij 


$  >  + 
o 


n^o 


$  >  <  $ 
n  n 


H’  -  VT.  ¥  > 
iJ  1  0 


E 


n 


(2.8) 


where  V?.  =  V(r..)  9  (£..  -  r..) 
ij  ij  iJ  ij 


s 


Define  the  operator  G^.  by 


G7.  Y  .  .> 
ij  1  iJ 


V,  .  h?  > 
ij  1  o 


(2.9) 


then  from  (2.6),  (2.8)  and  (2.9)  it  follows  that 


y  >  =  $  >  + 

o  1  o 


$  >  <  $  |  1/2^  ,  G?  . 

n _ n1  i,j  ij 


n^o 


E  -  E 


n 


$  > 
o 


+  .  .  . 


(2.10) 


Then  with  the  approximation  %  |$Q>,  we 


have  the  integral  equation  for  G 


ij 


<$  I  GT,  I  <f>  >  =  <  $  V,  .  $ 

o  1  ij  1  o  o  1  ij  1  o 


>  +  Z 

n^o 


<  $  I v? . I  $  >  <  $  |g?.|$  > 

0  1  ij  1  n  n 1  ij  1  0 


E  -  E 


n 


;  y  , 
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If  now  the  separation  distance  is  specified  by 
imposing  the  condition 


<$  | GT . I  $  >  =  0. 

o  '  ij  1  o  J 


(2.12) 


then  with  the  approximation  >  %  l$0>>  the  energy  shift 


is 


<$  H»  y  >  =  -i—  2_  <$  V? .  y 
o'  1  o  2  ..  o 1  ij  1  o 


1 


o  H  <$  | G®. | 

2  ij  u1  u 


■i-H  <*  |o!J*  > 

2  .  .  o  i  1  o 

ij  J 


=  os 

and  the  total  energy  of  the  system  is 


e  =  e  =ZL  <*  It.  +  1/2  21  vh  I*  > 

0  i  01  j/i  1J  0 


=  2.  E.  -  1/2  I.  <#  |V£ 
i  1  ij  °  1 


$  > 
0 


£  -> 

where  V.  .  =  V(r..)  6  (r..  -  £..)  is  the  long-range  part 
ij  ij  ij  ij 

of  the  interaction  potential  and  is  the  single-particle 


energy 


,  .  o  . 

E.  =  <3>  T .  $  >  +  Z_  <  $  V .  .  $  > 

1  0 1  1 1  0  jyj  o 1  lj  1  o 


Equation  (2.10)  can  be  compared  with  the  complete  reaction 
matrix  which  in  the  notation  used  by  Moszkowski,  Scott  has 
the  following  form 

G  =  V  +  V  G 
e 


16 . 


G | Kq>  (2.13) 

(2.14) 

Here  SQ,  t  denote  the  initial, intermediate ,  and  final 

relative  momentum  of  the  two  interacting  nucleons,  and  p 

is  their  average  momentum  so  that  the  laboratory  momenta 

of  the  two  nucleons  in  the  intermediate  state  are  p  +  S’  and 

p  -  .  | ¥n  (S  ,  p)  >  is  the  wave  function  for  two  nucleons 

interacting  through  the  potential  V(r). 

For  brevity,  the  parameter  p  will  be  omitted 

g 

hereafter.  From  the  definition  of  G  ,  it  is  seen  that 
equation  (2.12)  may  be  written  as 

<  t  IVs  I Y(k  ) 
o'  1  o 


with  Q(p , K ’ )  =  0 


If  K2  +  p2/4  <K2 


=  1 


K ' 2  +  p2/4  -  K2f 


If  Kr  -  p  / 2  >Kf 
—  otherwise 


K'p 


The  matrix  element  of  G  is 


<  K  G  K  >  =  <  K  V  |  K  >  +  <  K|  V  |  K T  >■ 
1  1  o  1  1  o  11 


Q 


e  ( K  ’  ,  K ) 


<K 


The  matrix  elements  of  G  are  also  given  by 


^  I  A  I  ^  -> 

<  K  G  K  ,  p  > 
1  1  o 5  ^ 


<  t|v|Tn(  £o>  p)  > 


> 


0. 


'  'r  i  di 

■  • 
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Thus  the  diagonal  elements  of  G  are  given  by 


<  t  | G |  t  >  =  <  1  | | Tn (S  )> 

o  1  1  o  o  1  1  o 

with  approximation  |  T >  =  | 'Fn  >  being  made.  It  is  well  known(l6) 
that  at  large  distances  (i.e.  r  >  £)  the  free  particle  wave 
function  is  a  good  approximation  to  the  two  particle  wave 
function  in  nuclear  matter. 

Thus  to  first  order  we  have 

<  t  | G |  k  >  »  <  t  |V£|  1  >  (2.15) 

o 1  1  o  o 1  1  o 


To  obtain  the  correct  integral  form  of  the  reaction 
matrix  we  have  to  take  into  account  the  spin  and  isospin 
states  of  the  two  nucleons,  which  are  required  for  the  com¬ 
plete  specification  of  the  dynamical  state  of  the  two- 
particle  system.  In  doing  so,  the  reaction  matrix  can  be 
written  as 


<  t  |G  I  t  > 


H  (2J  +  1) (2T  +  1) 
J ,  T 


r  oo 

J  K 


j  t(Kr)V®’ J(r )r2dr 

(2.16) 


where  is  the  potential  depending  on  the  total 

Ji>  J 

momentum  J,  orbital  angular  momentum  £,  spin  S  and 
T  of  the  interacting  nucleons. 


The  total  energy  of  the. system  is  then 


angular 

isospin 


E  =  It.  +  1/2  SI  <i?|G|t> 
i  1  1  j  J 

K.2 

=  £  +  1/2  SL.  <  t  |G  |  t  > 

i  i  j  j 

=  £  [  E (L)  -  1/2  U(t.)] 


(2.17) 
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where 


U(ti) 


K  =  1/2  (K.  -  K. ) 

i  3 

<  S  |  G  |  t> 

K.<Kf 

The  reaction  matrix  <  j£  |  G  |  t  >  is  defined  in 

terms  of  the  relative  momentum  It  of  the  particles  i  and  j  . 

To  find  the  average  potential  felt  by  the  particles  with 

momentum  ,  we  calculate  the  average  value  of  G  when  is 

kept  fixed  and  K.  takes  on  all  the  values  in  the  Fermi  sea. 

J 

Then  we  mutiply  the  result  by  the  average  density,  p=  2Kj; 


to  get 


U(Kj_) 


=  P 


<K I G I K  >  d JK  . 


Kj  <K; 


d3K 


KJ  <Kf 


3" 


(2 


The  average  potential  energy  of  the  particle  in  the  ground 
state,  U,  can  be  written  in  terms  of  U(K^)as 


U  = 


U(K. )dJK. 

'K1<Kf  1  7  J  Ki<Kf 


d3K 


or  directly  as  an  integral  of  G 
U  = 


PK 

=  p  /  1  dK  p(K,Kf)  <  K  |  G  |  K  > 


.18) 


2.19) 


(2.20) 


.  ■ 
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where 

P  (K,Kf) 


24  K2 


(1  - 


3K 

2Kf 


+ 


K 


3 


2K 


3 


) 


is  the  normalized  probility  distribution  of  relative 
momentum  for  two  particles  in  the  Fermi  sea.  For  low 
densities,  instead  of  equation  (2 . 20)  for  U,  we  can  use 
the  reaction  matrix  evaluated  at  the  average  relative 
momentum  K  multiplied  by  the  density.  The  average  relative 
momentum  is  a  function  of  which  depends  on  the  form 
of  <  K | G |  K> .  Usually  one  uses  K  =  yo . 3  which  is  exact 
if  <K  | G |  K>  is  a  simple  quadratic  in  K.  Then  U  becomes 


U  =  — Kf3<  K  |  G  |  K  > 
3  IT 


(2.21) 


The  average  kinetic  energy  T  of  the  particle  in 
the  Fermi  sea  is 

T  =  (3/10)Kf2.  (2.22) 

Thus  the  binding  energy  W  per  nucleon  in  nuclear 
matter  in  the  first  order  is  given  by 

W(Kf)  =  3/10  Kf2  +  Kf3  <  K | G | K  > 

3  u 


(2.23) 
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or  in  integral  form 


W(Kf)  =  3/10  Kf2  +  — L 

6  IT 


3 


S.  ( 2J+1 ) ( 2T+1 ) 
JST 


r  oo 


.2  /  tt  \  ttS  ,  T  /  \  2  , 

J £  (Kr)V£5 j(r)r  dr 


(2.24) 


The  integrations  of  the  last  part  of  (2.24)  are 

easy  to  carry  out  for  £  <_  2 .  For  £  >2  the  contribution 

( 2 ) 

to  the  reaction  matrix  is  very  small  so  we  believe  that 
it  is  good  approximation  to  neglect  the  contribution  for 
£  >2. 


The  solution  of  the  reaction  matrix  is  now  reduced 

to  the  problem  of  calculating  the  separation  distance  £ 

in  a  self-consistent  way.  A  first  order  calculation  can 

be  based  on  the  knowledge  of  the  two-nucleon  scattering 

phase  shifts  and  the  long-range  part  of  the  potential. 

The  first  order  approximation,  £  ,  is  obtained  by  letting 

the  interaction  energy  denominator  tend  to  zero.  Also 

the  centre  of  mass  of  the  two  interacting  particles  is 

assumed  to  be  at  rest  and  the  exclusion  operator,  Q,  is 

s 

replaced  by  the  principle  value  operator.  Then  Gp  ,  the 

g 

first  order  approximation  to  G  satisfies  the  equation 

Gp  =  Vs  +  Vs  -Ijr  Gp 
P  2  P 

where  e  ( K  , K 1 )  =  K  -  K1  is  the  first  approximation  of 
o  o 

eN  and  P  stand  for  Cauchy's  principle  value.  This  is 

just  the  equation  for  the  reaction  matrix  with  the  potential 
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g 

V  and  the  solution  Is  given  in  terms  of  the  scattering 
phase  shifts  due  to  the  short-range  potential  by 

<  K  |  G®  |  K  >  = - ^ -  2  (21  +  1)  C.  S.(K,  5  )  (2.27) 

K  i 

where  C  =  3/4  for  even  i  and  5/4  for  odd  i. 

Xy 

Prom  the  theory  of  scattering  it  is  expected  that 
the  limit  of  the  reaction  matrix  as  the  e^  -*e^  tends 

X/ 

g 

to  the  t-matrix  of  scattering  theory  (for  potential  V  ) 

which  is  given  in  terms  of  tan  5^  rather  than  . 

( 23  ) 

However,  Fukuda  and  Newton  have  shown  that  as 

energy  denominator  tends  to  zero,  the  limit  of  G?  depends 
on  as  shown  in  (2.27). 

In  the  case  of  spin-  and  isospin-dependent 

g 

potentials  the  Gp  matrix  is 

<  K  I  G®  I  K  >  =  Z  (  2  J  +  1 )  ( 2T  +  1)  6  ^ 5  T  (2.28) 

2K  ^  J 

O  m 

where  6 „ 3  ^  is  the  phase  shift  due  to  the  short-range 

S  T 

part  of  the  potential  V^5j  (r)  0  (  5Q-r).  We  require 
the  first  order  separation  distance,  £  (K),  be  specified 
by  the  condition  that  the  diagonal  elements  of  Gp  vanish 

<  K|g!;|K  >  =  0.  (2.29) 

i  p  i 

To  find  the  separation  distance  we  define  the 

ST  (24) 

phase  function  (K,r)  by  the  differential  equation" ' 


' 

■ 
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with  the  boundary  condition 


(2.31) 


where  6^5j  (K)  are  the  known  scattering  phase  shifts 


S,T 


S  T 

produced  by  the  complete  potential  V  5T  (r). 

36  *  J 


For  practical  purposes  one  prefers  to  integrate 


S ,  T 


6  ’  (K,r)  from  the  origin  to  the  point  £  which  is  de- 

X/  j  u  O 

fined  by  eq.  (2.29).  In  the  case  of  a  hard  core  potential 
of  radius  c,  one  starts  at  the  core  radius  with  the  condition 
that  the  wave  function  should  be  zero  inside  the  core.  This  im¬ 
plies  that  the  phase  function  at  the  edge  of  the  core  is 
given  by 


(2.32) 


tan  6^j  (K, c ) =  (Kc)/n&(Kc) 


One  can  also  integrate  the  first  few  partial  phase  function 

(e.g.  £=0,1  and  2)  with  no  phase  scattering  from  infinity  to 

point  £  where 
o 


(2.33) 


The  solution  £  (K)  can  be  used  to  compute  approximate 
values  for  the  diagonal  elements  of  the  reaction  matrix  and 
the  single-particle  potential  energy  from  eq.  (2.l6)  and 
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(2.18).  At  first  sight  it  may  seem  that  the  result  is 
nearly  the  same  as  calculating  Gp ,  i.e. 

<  K  |  Gp  |  K>  =  -  £1  (  2J  +  1 )  (  2T+1 )  (2.34) 

J  ,T 

But  this  is  not  the  case  and  the  difference  between  G 
calculated  from  (2.16)  and  (2.28)  may  be  large.  For 
example,  the  one-pion  exchange  (OPEP)  contributes  nothing 
to  the  G  defined  by  (2.16),  but  it  does  contribute  to  Gp 
in  (2.28).  It  should  be  noted  that  the  approximate  re¬ 
action  matrix  given  by  eq.  (2.16)  is  a  much  better 
approximation  for  the  actual  nuclear  matrix  than  Gp . 

The  reason  is  that  in  nuclear  matter  the  attractive 
potential  cannot  change  the  Fermi  gas  nearly  so  much  as  it 

(5) 

affects  free  particles.  Moszkowski  and  Scott  have  shown 

that  Gp  is  a  useful  approximation  only  for  large  K. 

To  find  the  separation  distance  by  eq.  (2.33)  one 

encounters  difficulty  for  the  coupled  states  due  to  the 

coupling  by  tensor  force.  For  weak  coupled  states  the 

tensor  force  can  be  neglected  but  for  the  strong  coupled 

states  such  as  S1  -  we  can  only  find  the  separation 

distance  by  solving  the  Schrddinger  equations. 

2  3 

The  Schrddinger  equation  for  JS1  -  'JD1  states  are 

,2  0 

H?  +  {K  -  vn(r))u  =  V12W  (2.35) 

dr 

+  (K2  -  -§-  -  v  (r)w  =  V12u  (2.36) 

dr  r 


• 
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where  u  is  the  wave  function  of  state  and  w  of  state. 

V.^2  is  the  tensor  force,  and  is  the  central  potential. 

3 

For^D^,  V^2  is  a  combination  of  central  *V  ,  spin-orbit, 

VLg } andtensor  potential,  .  The  explicit  forms  of  these 
potential  are: 


V 


11 


V 

c 


V12  =  2/2  VT 


3V 


LS 


The  Schrbdinger  equation  of  free  particle  in  S 
state  is 

2 

- - %  +  K2  =  0  (2.37) 

d  r 


where  <j)  is  related  to  the  spherical  Bessel  function  as  0  = 
rj  (Kr).  Equation  (2.37)  can  be  written  as: 


- — |  +  (K2  -  Vll(r»  * 
dr 


-VX1  (r)* 


Substracting  (2.35)  from  (2.38) 


(2.38) 


dr 


‘U  (4,  -  u)  +  (K2  -  Vl;L(r))(*  -  u)  =  (-V12w  -  V^)  (2.39) 


which  can  be  written  as 


^5-  X  +  (K2  -  Vn)  X  =  (-V-,  5w  -  V„*) 


dr' 


12 


11 


(2.40) 


where  x  =  <|>  -  u  is  the  difference  between  the  free  particle 
wave  function  and  the  perturbed  wave  function. 
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If  we  put  w  =  —  £  where  w  is  the  perturbed  wave 
function  of  state,  then  we  have  the  coupled  differential 

equation : 


- - b  +  (K2  -  V  )  X  =  (V i p  c 

d  r  11 


(K2  -  ^2  -  V22)  C  =  -V12  (♦  -X) 


(2.41) 


(2.42) 


with  the  boundary  condition 

X  (c)  =  4>  (c)  (2.43) 

C(c)  =  0  (2.44) 

The  first  order  separation  distance  is  a  point 
where  the  perturbed  wave  function  joins  smoothly  to  the 
unperturbed  wave  function.  Since  x  and  x’  are  continuous 
functions  of  r,  it  follows  that 

X  (e0)  =  0  (2.4) 

and 

X'(E0)  =  0  (2.46) 


The  numerical  results  are  given  in  Chapter  4. 


r  •  iv-.n  :  • 
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Chapter  3* 


First  and  Second  Order  Correlations 


Once  the  approximate  separation  distance  is 


obtained  from  (2.33),  we  can  determine  single-particle 
energies  which  enable  us  to  find  a  more  accurate  value  of 
the  separation  distance,  £,  by  solving  (2.10),  subject 
to  the  condition  that  <  V  |GS|  *F  >  =  0.  This  can  be 
done  most  easily  by  using  the  coordinate  representation 
of  the  wave  function.  Assuming  that  E  -  En  in  (2.7) 
depends  only  on  the  excitation  energy  of  the  two  particles, 
we  can  reduce  the  equation  to  one  describing  the  perturbed 
wave  function  of  the  two  interacting  particles: 


rather  weak,  and  also  because  of  the  exclusion  principle, 


the  distortion  of  the  eigenfunctions  of  H  ,  I  4» 0  >  ~e 


are  small  .  To  a  good  approximation  they  can  be  regarded 
as  plane  waves.  Thus 


'F  (t‘r)  =  exp  (iK*r)  +  - o 

(2ir)J 


I 


J 


KT>K 


f 


exp  IK  *  •  (r  ~r  )VS  (r  *  )  ¥  (  Kt  ) 

e  ( K T  , K ,p  ) 

(3.2) 


if . 
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which  is  similar  to  the  Bethe-Goldstone  equation. 

To  satisfy  the  condition  <Y  | Gs | ^  >  =  0,  the 

solution  of  (3*2)  must  satisfy  the  condition 

\exp(-ii?-r)Vsf (]?•?)  d3r=  o  (3.3) 

The  energy  denominator  in  (3.2)  is  given  by 
e(K,K',p)  =  E^.tj)  -  E(£mJn) 

=  E(|  p  +  t,  |  p  -  t)  -  E(|  p  +  K',  |  p  -  t') 

where  E(K^,Kj)  and  E(Km,Kn)  are  the  eigenvalues  of  the 
two-particle  Hamiltonian  Hq  defined  by  relation  similar  to 
(2.1)  and  (2.2). 

From  (3*2)  and  (3-3)  a  better  approximation  for  £ 
can  be  found  which  in  turn  will  determine  e(K’,K,p)  more 
accurately.  This  process  can  be  continued  until  self- 
consistency  is  achieved,  i.e.  until  the  difference  be¬ 
tween  two  successive  values  of  £  is  negligible.  Practically 
this  is  a  complicated  method  and  is  difficult  to  isolate 
and  determine  the  roles  of  the  physical  aspects  of  the 
problem,  such  as  the  effects  of  the  exclusion  principle 
and  the  motion  of  the  centre  of  mass. 

To  simplify  the  mathematics  we  introduce  the  method  of  re- 
ference  spectrum  of  Bethe  et  al ^  .  Eq .  (3-2)  and  (3*3) 

are  solved  under  the  following  conditions-  (a)  The  region 
of  integration  over  K’  is  extended  to  include  the  Fermi 
surface;  (b)  we  assume  that  the  centre  of  mass  of  the  two 
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interacting  particles  is  at  rest,  i.e.  p  =  0;  (c)  we 
fit  e(K,K’)  with  a  quadratic  form  as  a  function  of  Kf, 
e(K,K')  =  [K  +  y  (K)]/m#  where  m*  is  the  effective  mass. 

In  (c)  we  have  assumed  that  e  remains  a  negative  quantity 
for  all  values  of  K  and  Kf  and  this  implies  that  no  real 
scattering  can  take  place  in  Fermi  sea.  The  resulting 
partial  differential  equation  is 

2  2 

(V  -  y  )[y(J)  -  exp  (iS*r)]  =  m*V(r)  d(£-r)y  (y)  (3,.  4) 


In  the  case  of  a  central  potential  we  can  use  the 
partial  wave  decomposition  of  y  (?) 

oo 

4/  (J)  =  _i—  n  (2  £+l)i£P  (cos  e)u  (r)  (3.5) 

r  *  =  o  36  36 

to  reduce  equation  ( 3  •  *0  to  a  set  of  total  differential 
equations : 

[_1_  _  i(z.  +1)  _  y2][u  (r)  -  rj  (Kr)]  =  m*  Vs(r)u  (r)  (3.6) 

dr  r 

with  the  boundary  conditions  u^(c)  =  0  and  (3.7) 


l~d(uz  (r)  -  0  (Kp)  )/dp  -j  =  j-d  (rh& ( 1 )  ( jyr )  )/dr  j 

u  (r)  -  rj  ^ ( Kr )  rh£(iyr) 

r=C  r 

where  h  is  the  spherical  Hankel  function  of  the  first 

X/ 

kind . 

The  separation  distance  is  defined  by  equation  (3.3) 


(3.8) 


which  can  be  expressed  in  partial  waves  as  follows: 


. 


' 


.  1 .  s 


■ 
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TT 


o-  (2J+1)(2T+1) 
J,T 


'C(K) 


rj£(Kr)V(r)  u£(r)  dr  =  0 


Substituting  the  integrand  for  V(r)u  (r)  from  (3.6),  we 
find 


TT 


2 


2_  ( 2J+1 ) ( 2T+1 ) 

J  5  T 


■C(K) 


rj£(Kr)V(r)  u£(r)  dr 


TT 


>_  ( 2J+1 ) ( 2T+1 ) 

J  3  T 


1 


m 


* 


.  /t,  s  ,  d2  1(1  +  1) 

rU(Kr)(-— 2  -  - - 2-1 

^  o  dr  r 


__  *> 


( u  ^  -  rj£(Kr))dr 


2_  (2J  +  1)  (2T+1)  |(:y2  +  K2)  |  rj£(Kr)  (rj£(Kr)  -  u£(r)j 

J  5  T 


+ 


[rj£(Kr) 


d 


■{  u  (r )  -  r j  ( Kr ) }  -  {  u0(r)  -  rj.(Kr)} 


dr  L  L 


SF-  rh(Kr) 


=  0. 


r  =  C 


The  last  step  involves  integration  by  parts  and  substitution 
of  the  boundary  condition  (3.8).  As  mentioned  before,  at 
the  region  r  >  £  it  is  good  approximation  to  replace  the 
exact  wave  function  by  the  plane  wave.  In  doing  so  the  last 
term  of  eq .  (3*10)  can  be  neglected. 


(3.9) 


dr 


(3.10) 
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Then  (3.10)  becomes 


*  f%(  K) 

— p—  Z~  (2J  +  1)(2T+1)  /  rj  (Kr)  V  (r)  u  (r)  dr 

J ,  T  Jo  1 

7T  "  p  P  /"£  (  K  ) 

=  -  X.  (2J+1) (2T+1) (y^+Kj  / 

2m*  J  i  T  JQ 


rj  £(Kr )x£(Kr)dr 


(3.11) 


where  X£(Kr)  =  rj£(Kr)  -  u£(r). 

Using  the  approximation  u  (r)  ^  rj  (Kr)  for  r>  £  in  right 

X/  X/  0 

side  of  (3.11)  again  we  get 


TT 


(2J+1)(2T+1)[* 

J ,  T  m* 


2  2 


VK> 


J 


r  j  £  ( Kr )  x  £  ( Kr ) dr+ 


o 


(% 


r2j2  ( Kr )V(r )dr ] =0 . 


(3.12) 


For  hard  core  potential  the  exact  wave  function  inside  the 
core  is  zero,  i.e.  u£(r)  =  0  for  r  <c,  so  that  (3.12)  can 
be  written  as 


H  C  2J+1 )  (  2T+1 )  [j  r2j  2 

J,T  m  J„ 


(Kr)dr  + 


r 


2  . 

r  J 


£ ( Kr ) x  £ (r )dr ] 


C? 


H(2IT+1)(2T+1)  r2  j  p  ( Kr  )V(r  )dr  . 
2  J ,  T  ^ 


(3.14) 


The  quantity  y  (K)  is  related  to  the  single-particle 
potential  energy  u(K1),  which  can  be  written  as  quadratic 

2 


function  of  k. . 

l  ^ 

U(k.)  =  (-i-  -  1)  — 


+  A 


m 


* 


a+bK-j 
ki 


9 ( Kf  -  Kj) 


1+cK-i 2  2 


(3.15) 


. 
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where  A,  m* ,  a,  b  and  c  are  constants , and  0  is  the  step 
function,  then  the  energy  denominator  can  be  written  as: 


e ( K, K ’ ) 


and  thus 


Y2(K) 


2(2  ? 
K  -  K  ^  _  a  +  b  n 

m*  1  +  cK2 

m* ( a  +  bK2)  v2 

2  ^ 

1  +  cK 


(3.16) 

(3.17) 


Using  the  first  order  approximation  £  to 

o 

calculate  the  single-particle  potential  energy  u(K^) 
and  fitting  in  (3*15)a  (3*16)  and  (3.17)  we  can  have 
the  first  approximation  to  the  energy  denominator  e(K,K') 
and  the  effective  mass  m*  .  Then  we  calculate  x»  by 

X/ 

solving  eq.  (3.6)  and  integrate  (3.14)  to  find  a  more 

accurate  separation  distance  ?.  This  process  can  be 

continued  until  self-consistency  is  achieved. 

For  the  particle  above  the  Fermi  Sea 

the  contribution  of  odd  states  is  very  small  and  can  be 
(29) 

neglected  .  Thus  we  can  use  the  same  process  used  for 
particles  inside  the  Fermi  sea  to  calculate  single-particle 
energy  above  the  Fermi  sea  but  only  take  the  even  states 
with  statistical  weight  1  instead  of  3/4  into  account. 

After  calculating  the  first  order  separation  distance  and 
hence  the  first  order  potential  energy,  the  second  order 
approximation  can  be  calculated  easily. 

Once  the  exact  separation  distance  is  found,  the 
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contribution  of  the  long  range  potential  may  be  obtained 
from  the  Born  series.  Thus 


V 


£ 

exact 


=  V 


+  V 


Q 


V 


(3.18) 


If  instead  of  the  plane  wave  approximation  we  were  using 
the  unperturbed  wave  function  |  <j>  > ,  we  could  write 

.£  ^  ^ 


V"  -  <r  *K-r  iv* 

exact  e 


(3.19) 


where  the  integral  equation  for  |  <J>q>  is  given  by 


(Kr)  =  e 
o 


iK-r 


+  dr  ’  e 


it  •(?-?')  Q 


e  V  (r»)  ♦<)(Kr»)  (3. 


By  successive  iteration  of  (3.20)  and  substitution  in 

(3.19)  we  find  (3*18).  There  will  be  some  scattering  to 

states  outside  the  Fermi  sea.  The  contribution  from  this 

scattering  is  expected  to  be  quite  small  for  the  following 

£ 

reasons:  (1)  V  is  already  small,  and  (2)  the  Pauli 

£ 

principle  strongly  inhibits  any  scattering  caused  by  V  . 

£  Q  £ 

The  diagonal  elements  of  V  — —  V  for  each  state  can  be 
written  explicitly  as 


K> 


S ,  T 
£ ,  J 


(■»  <K  |  vV’>  <K’  |V£|K>  S’J  Q(p,K 

Jo  e(K,K',p) 


It ii  K'2 

(  2ir )  3 


dK' 


20) 


I 


■ 


. 
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Q(p.K' ) 


,2 


o  e  ( K,  K  '  ,  p  )  (2  7T ) 


1,71  K'—  d  K'[l  j£  ( K '  ,r)V(r)2’J  j ^ ( Kr )r 2dr ] 2 


=  8 


o  e 


Q'f >K'  }  K'2dK’  [f  j  (K'r)V(r)®’J  j  (Kr)r2dr] 
(K,K'p)  Jc  1  l’J  l 


(3.21) 


then  for  all  partial  waves  this  contribution  is 


<K|V£  V£  |K>  =  -g-  L  ( 2J  +  1)  ( 2T+1  )<K  |  V£  V£|K>S,T 


j3t 


=  >_  ( 2J+1 ) ( 2T+1 ) 
J,T 


C  oo 


^  v  K,2dKT 
e  (K  j  K ' p  ) 


l  j  J 
(3.23) 


j  £ (K ' r )  V£ 5  J  j£(Kr)r2dr]2 

where  p  is  the  centre-of-mass  momentum  of  the  interacting 
pair . 


The  effect  of  p  is  greatest  in  the  exclusion 
operator  Q  and  is  not  so  important  in  the  energy  denominator. 
For  this  reason  we  replace  p  in  Q  by  the  average  (over  the 
Fermi  sea)  centre  -  of  -  mass  momentum  of  two  particles 
having  a  relative  momentum  K.  Thus 


Kf(Kf 


K) 


1+  - - - 1  for  K  <  K 

L  3Kf(2Kf+K)  -J 


(3.23) 


=0.  K  >Kf 

A  cruder  approximation  was  used  in  the  energy  denominator 

where  the  value  p  =  was  used.  This  value  corresponds 

to  the  peak  in  the  distribution  of  the  number  of  pairs 

( 5 ) 

of  particles  having  centre  -  of  -  mass  momentum  pw  . 
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Chapter  4. 


Numerical  Calculation  with  BGT  Potential 


A.  Details  of  Calculation  and  Results 

The  Method  described  in  the  previous  chapters  has 
been  used  to  calculate  the  effective  interaction  for  the 
Brueckner  -  Gammel  -Thaler  (BGT)  potential.  The  BGT  potential 
contains  a  central,  a  spin  -  orbit  and  a  tensor  term  and  has 
Yukawa  well  outside  a  hard  core  . 


V(r)  =  + 
V(r)  =  V 


e  yTr 


-y  r  -u, nr 

e  c  f  17-  e  LS  ^ 

c  ucr  •  LS  pLSr  §12VT 


The  parameters  of  this  potential  are  shown  in  Table  I. 


TABLE  I. 


r<c 

r>c 

(4.1) 


Parameters  of  the  BGT  Potential 
The  Core  Radius  is  0.4  Fermi 


SPIN 

PARITY 

V(MeV) 

IlfF'1) 

Central 

Singlet 

even 

-434 

1.45 

singlet 

odd 

+130 

1.00 

triplet 

even 

-877. 4 

2.091 

triplet 

odd 

-14 .0 

1.00 

Spin-Orbit 

Triplet 

even 

-5000 

3.70 

triplet 

odd 

-7315 

3.70 

Tensor 

triplet 

even 

-159.4 

1.045 

triplet 

odd 

+  22.0 

0.08 

•IV 

i 
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TABLE  II, 


Phase  Shifts  Calculated  from  the  Phase  Functions 


The  Phase  'Shifts  are-in  radians 


E(MeV) 

!s 

^0 

1 — 1 

tu 

1 — 1 

OJ 

Pi 

1 — 1 

3p 

0 

3P  3 
1 

P  ^ 

2 

D2 

1.1 

30 

.789 

-.14 

.0154 

.248 

-.112 

.08 

.0764 

1.2 

35.8 

.741 

-  .16 

.021 

.  26 

-  .12 

.09 

.10 

1.3 

42 

.693 

-.18 

.027 

.269 

-.14 

.11 

.128 

1.4 

48.7 

.646 

-.20 

.0358 

.271 

-.15 

.13 

.158 

1.5 

55.9 

.6 

-.218 

.0446 

.27 

-.17 

.148 

.19 

1.6 

63.6 

.55 

-.236 

.054 

.26 

-.18 

.167 

.224 

1.7 

71.9 

.511 

-.25 

.0649 

.255 

-.198 

.185 

.259 

. 
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TABLE  III. 


Phase  Shifts  Calculated  by  Gammel  &  Thaler  from  4205  Potential 

The  Phase  shifts  are  in  radians 


(MeV) 

L 

o 

% 

3p 

o 

rH 

Ph 

OO 

3? 

2 

30 

.771 

.014 

.219 

-.121 

.024 

40 

.695 

.  024 

.2365 

-.149 

.0995 

60 

.  5068 

.  047 

.2365 

-.194 

.14555 

80 

.466 

.0725 

.213 

-.231 

.1835 
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Pig. Ia  -  The  Phase  Function  of  State  as  a 
Function  of  r(F)  for  =  1.5  F 
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Phase  Functions  (radians) 


40. 


Fig.  Id  -  The  Phase  Functions  of  ^P^  and  ^P^ 
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The  spin-orbit  and  tensor  operators  £•§  and  §^2  have 
their  standard  meaning,  i.e. 

£•3  =  [J(J+1)  -  L(L+1)  -  S(S+1)]  (4.2) 


L’L” 

J-l 

J 

J+l 

J-l 

-2( J-l ) 

0 

6[J( J+l ) ] 1/2 

2  J  +  l 

2J  +  1 

VT(r)  -J 

0 

2 

0 

J+l 

l h 

6CJCJ+1)] 

0 

-2 ( J+2 ) 

2J  +  1 

2J+1 

This  matrix  relation  can  be  obtained  by  acting  the 
tensor  operator  $^2  on  the  angular  and  spin  dependent  wave 
functions . 

Before  finding  the  first  approximation  to  the 
separation  distance, 5  ,  we  integrate  (2.30)  from  the  core 
radius  with  the  boundary  condition  (2.32)  to  infinity  to 
get  the  phase  shift, 6  *  (K),  the  result  is  shown  in 

A/  ^  cJ 

Table  II.  The  result  has  been  compared  with  those  calculated 
by  Gammel  and  Thaler^2^  using  GT4205  potential  (See  Table 
III),  and  shows  the  excellent  agreement  with  the  direct 
integration  of  Schrddinger  equation.  Figs.  1  are  the  phase 
functions  for  1SQ , 1P1 , 3Pq ,3p1 , 3p2 , 1d2 ,  and  3D2  states  as 
functions  of  r(F)  . 

In  all  of  our  calculations  we  used  the  fourth-order 
Runge-Kutta  method  to  solve  the  first  order  differential 
equation  and  the  Simpson’s  rule  for  integration.  Due  to 
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the  strong  coupling  by  the  tensor  force  we  calculate 

the  separation  distance  for  state  by  solving  the 

coupled  differential  equation  (2.35)  and  (2.36)  as  stated 

in  Chapter  3.  The  coupling  between  and  is  small, 

so  that  we  neglect  the  coupling  in  state.  We  also  cal- 

culate  the  state  separately  because  its  separation 

distance  is  close  to  the  core  radius  and  the  contribution 

of  to  the  potential  energy  is  small.  Then  it  is 

easy  to  integrate  the  phase  function  (2.30)  for  all  the 

33  3 

states  £<2  except  S^,  and  starting  from  the  core 

radius  to  the  point  where  equation  (2.32)  is  satisfied. 
Moszkowski  et_  ai^ave  pointed  out  that  the  separation 
distance  depends  upon  the  momentum  and  also  on  the  angular 
momentum.  However,  in  our  calculation,  we  introduce  only 
one  separation  distance  for  all  partial  waves  (except  those 
mentioned  above).  The  result  of  our  calculation  shows  that 
the  approximate  separation  distance  is  nearly  independent  of 
and  K,  and  it  has  a  value  of  1 . 07F( 1F=10-I^cm) .  This 
result  (i.e.  £  independent  of  K)  may  be  related  to  the  fact 
that  the  cross-section  for  nucleon-nucleon  scattering  re¬ 
mains  constant  over  a  wide  range  of  energies.  After 
finding  the  separation  distance,  the  potential  energy  can 
be  easily  calculated  from  (2.16)  and  (2.21)  by  using 
numerical  integration. 
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3 

To  calculate  the  separation  distance  for  state, 

we  factorize  the  second  order  differential  (2.4l)and  (2.42) 

( 1 ) 

into  the  following  six  first  order  differential  equation  . 


1 

=  -  [  1  +  f(r)  Y.^] 

(4.4) 

dr 

dY2 

=  -  Y1[f(r)  Y2  +  Vn(t,  -  V12  C] 

(4.5) 

dr 

dY3 

=  f(r)[Y3  Y1  +  Y2] 

+tVll*  -  V12  ^ 

(4.6) 

dr 

dZ1 

=  -[1  +  g(r )  Zx2] 

(4.7) 

dr 

dZ  2 

=  -Zl[v12(*  -  x)  + 

g  (r)  Z2] 

(4.8) 

dr 

dZ^ 

=  g  (r)[Z3  Zx  +  Z2] 

+  v12  U  -X  ) 

(4.9) 

dr 

Here  f(r) 

’  k2  "  V11 

(4.10) 

g(r ) 

n  2  22 

r 

(4.11) 

4>  ( r ) 

=  r  j  0  ( Kr ) 

(4.12) 

X 

=  Y  Y  +  Y 
*1  X3  2 

(4.13) 

And  ; 

=  Z1  z3  +  Z2 

(4.14) 

Equation  (4 

.  4)  and  (4.5)  are  subject  to  the  following 

boundary  conditions: 

Yx(c)  =0 

(4.15) 

Y2(c)  =* 

(c) 

(4.16) 

^5. 


The  point  where  Y^Cr)  is  zero  is  the  desired  point  , 

i .  e . 

V  eo)  =  0 


At  this  point  we  put 


Y3(c0)  =  0 


(4.17) 


(4.18) 


Since  x  and  x’  are  continuous  functions  of  r,  it  follows 
that 

X  UQ)  =  x’  (50)  =  0 

Assuming  the  continuity  of  5  at  r  =  e,  ,  we  have 

0 


(4.19) 


S'  (  Sn) 

4  <  ?0) 


d[rj  ,?(Kr)] 


dr 


r j  ^ ( Kr ) 


(4.20) 


r  =  C, 


To  start  integration  we  use  the  trial  function^  =  -rJ2(Kr)  and 
integrate  (4.4)  and  (4.5)  from  the  core  radius  with  the 
boundary  conditions  (4.15)  and  (4.16)  outward  to  the  point 
E,  satisfying  (4.17),  then  with  the  boundary  condition  (4.18) 
we  integrate  (4.6)  inward  to  c  and  calculate  x  from  (4.13). 
Knowing  x  we  are  able  to  integrate  the  equation  for  and 


from  c  to  £  with  the  boundary  conditions 

Z1(c)  =  0 
Z2(c)  =  0 

From  equation  (4.7),  (4.8),  (4.9),  (4.14)  and 


(4.211 

(4.22) 


(4.19)  it  follows  that 


Z3(5o}  “  - 


(  )  (  d^J?(KrU 

2  0  V  dr 


r  =So- 


(4.23) 


r  j  p  ( Kr )  +  Z  (e  ) 

d  0  dr 


r  =  ? 
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With  this  value  of  Zq  we  can  integrate  (4.8)  from  E, 

j  o 

to  c  and  find  c  from  (4.14).  Usually  two  or  three  iter¬ 
ations  are  enough  to  give  an  accurate  value  of  £  ,  y 

o 

and  c  • 

The  separation  distances  for  state  are  given 
in  Table  IV.  It  is  interesting  to  note  that  the  separation 
distance  for  does  not  depend  strongly  on  K  at  low 

energies . 

The  potential  energy  for  state  can  be 

calculated  as  for  'v  states.  The  separation  distance 
3  3 

for  and  are  very  close  to  core  radius,  thus  we  use 

the  following  approximation  to  calculate  the  potential 
energies:  3 

u,  =  - 4 -  •  —  ( 2J+1 ) 6  J  (K,  °0  (i).24) 

^Dj  3-rr^  2  1 

The  potential  energies  for  all  of  the  states  are 
given  in  Table  V. 

TABLE  IV. 

The  Separation  Distances1  for  S-.  State 


Kf(F  -1)  1.1  1.2  1.3  1.4  1.5  1.6  1.7 

p  (F)  .97  .98  .99  1.01  1.02  1.05  1.09 

0 
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To  calculate  first  and  second  order  corrections 
we  have  to  calculate  also  the  potential  energy  above  the 
Fermi  sea  and  fit  the  potential  energy  in  and  above  Fermi 
sea  into  equation  (3.15).  Then  from  (3-15),  (3.16)  and 

(3*17)  we  can  get  the  values  of  the  effective  mass,  energy 

2 

denominator  and  7  (K).  The  single  particle  potential 
energy  as  function  of  relative  momentum  for  =  1.5F~^ 
is  shown  in  Fig  2.  The  values  of  the  constant  for  equation 
(3*15)  are  shown  in  Table  VI. 


TABLE  VI. 

2 

The  Values  of  Effective  Mass,  y  (K)  and  the  Constants 

for  Energy  Denominator 


Kf(F  1) 

2  -? 

Y  (F  L 

m* 

<AJ 

1 

Pc* 

b 

c(F2) 

1.4 

1.074 

.95621 

3.25947 

-.01132 

1.5476 

1.5 

1.1739 

.94666 

4.00285 

-.01393 

1.5476 

1 . 6 

1.4206 

.  93603 

4.85798 

-.0169 

1.5476 

m*  and  b 

are  dimensionless 

numbers . 
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K(F 
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The  first  order  corrections  to  the  separation 
distance  and  potential  energy  are  calculated  from  (3.6) 
and  (3ol4)0  To  calculate  x  we  solve  the  second  order 

X/ 

differential  equation  by  the  method  similar  to  that  used 

o 

in  the  calculation  of  £  for  JS, .  We  factorize  (3.6) 

o  1 

into  three  first  order  differential  equations  as  (4.4), 
(4.5)  and  (4.6),  but  replace  f(r)  by 


_  (  +  y2  +  m*  V) 


and 

V11  cj)  -  V-,2  t  by  -m#  V  r  j  %  (Kr). 


The  boundary  conditions  of  (2.45)  and  (2.46)  are  also 
replaced  by 


X 


(  5o) 


d[rh^  (iy  r)]/dr 
rh^  (iy  r) 


(4.25) 


r  =K 


and 


{d[rh^  (iy  r)]/dr)  Y 2 

{d[rh^'j  (iy  r)]/dr}  Y^rh^  (iy  r) 

(4.26) 


where  (4.26)  is  obtained  from  (4.25). 

The  first  order  correction  of  separation  distances 

and  contribution  potential  energies  are  given  in  Table  VII. 
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The  second  order  potential  terms  can  be  obtained 
by  integrating  (3.22)  with  corrected  separation  distance  £. 
The  result  is  given  in  Table  VIII. 

TABLE  VIII. 

The  Contribution  of  the  Second  Order  Terms  to  the 

Potential  Energy.  All  unmarked  Units  are  in  MeV. 


Kf(F  L 

h 

0 

i — 1 

OQ 

on 

i — 1 

P-I 

i — 1 

3P 

0 

1 — 1 

P-. 

on 

3P 

2 

\ 

3°2 

TOTAL 

1.4 

-.143 

-.05 

-.267 

-.024 

-  .  656 

-.689  -.295 

00 

-=3“ 

1 

-2.61 

1.5 

-.2 

-.07 

-.209 

-.033 

-.61 

-.684 

-.34 

-.525 

-2.67 

1 . 6 

-  .454 

-.14 

-.151 

1 

• 

0 

-.527 

-.635  -.375 

-.544 

CO 

00 

C\J 

1 

The  binding  energies  per  particle  in  nuclear  matter 


is  given  in  Table  IX. 
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B.  Discussion. 

(2) 

The  calculations  of  Brueckner  and  Mastersonv  '(BM) 
and  of  Razavy  both  give  about  8  MeV  per  particle  binding  energy 
for  equilibrium  density  with  k^  =1.28  and  1.12f-^  respectively. 
G.E.  Brown  e_t  aJ  have  investigated  some  possible  causes  of 

this  discrepancy  very  thoroughly.  They  found  that  in  both 
calculations  the  binding  energy  was  overestimated,  making  the 
discrepancy  worse.  In  BM  calculation,  an  angle  average  was 
used  and  the  average  total  momentum  was  inserted  in  the 
exclusion  principle.  However,  this  is  not  a  good  approximation 
and  causes  an  error  of  about  1.7  MeV.  In  addition,  BM 
calculation  is  approximate  because  of  incorrect  treatment  of 
the  off-energy  effects  in  intermediate  states.  They  neglected 
the  effective  mass  introduced  by  off-energy  shell  propagation 

which  causes  an  error  of  2-4  MeV/particle .  Razavy,  on  the 
other  hand,  overestimated  the  binding  energy  by  omitting  the 
Pauli  and  spectral  corrections.  Correcting  for  this  would 
lower  the  binding  energy  by  about  the  same  amount  as  that  of 
Brueckner  and  Masterson, 

In  our  method,  by  adopting  the  intermediate  energy 
spectrum  from  the  work  of  Bethe  et  al.,  we  can  include 
higher-order  corrections  in  our  work  as  has  been  done  in  theirs. 

5 

Besides,  the  short-range  part  of  the  reaction  matrix  G  in  our 
work  is  of  Bethe  et  al . ,  with  the  exception  that  here  the 
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energy-gap  parameter  has  not  been  introduced  directly  and 

the  role  of  the  hard  core  in  determination  of  the  energy 

spectrum  has  been  reduced.  In  the  methods  of  Moszkowski  - 

Scott  and  Bethe,  the  short  range  reactions  matrix  makes  a 

positive  contribution  which  cancels  part  of  the  long  -  range 

q 

effect,  while  in  our  work  G  is  zero,  and  we  get  a  larger 
separation  distance  than  theirs,  hence  a  weaker  potential 
outside  and  a  better  convergence  for  the  perturbation 
expansion . 

The  result  of  our  calculation  shows  that,  instead 

of  Kf=1.5  F  \  the  saturation  point  is  at  K.^  1.43F  ^  * 

2Kf2  -1 

closer  to  the  correct  density p  =  — —  with  Kf=  1.4F 

3  it1*  1 

The  binding  energy  at  this  density  is  10  MeV/particle .  (see 
Fig.  3).  We  have  only  calculated  the  first  iteration  for  the 
separation  distance  .  Since  the  total  correction  to  the 
potential  energy  is  small  compared  with  the  first  order,  we 
believe  that  this  is  a  good  approximation.  The  discrepancy 
between  our  calculation  and  observed  value  is  still  about  5  MeV/ 
particle . 

Recently  much  work  has  been  done  to  account  for  this 
discrepancy,  Bethe^^)  has  studied  higher  order  diagrams 
in  Goldstone  expansion.  His  estimate  in  third  order 
diagrams  shows  that  they  are  attractive  and  indeed  very 
important.  His  estimate  is  that  the  third  order  diagrams 
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Pig.  3.  The  Kinetic  Energy,  Potential  Energy 


57. 


approximately  cancel  the  correction  of  Brown,  Schappert, 
and  Wong^2^  and  bring  the  energy  back  to  the  value  of 
Brueckner  and  Master  son.  Moszkowski^^  have  also  found 
that  the  short  range  three  body  correlations  will  contribute 
to  binding  energy  less  than  2  MeV/particle . 

The  binding  energy  with  the  third  order 
corrections  is  still  smaller  than  the  observed  15.5 
MeV.  To  obtain  a  larger  binding  energy,  Wong^2^  has 
suggested  to  replace  hard  core  by  a  soft  repulsive  core. 

For  Breit’s  potential  this  contributes  about  1.5  MeV. 


» 


- 
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Chapter  5. 

Variational  Method 


Consider  the  total  Hamiltonian  of  a  system 
with  A  particles  containing  the  sum  of  kinetic  energy  of 
the  particles  and  two-body  interactions  which  is  given 
by  _A_ 

h  =  n  t  +  -i-  v(?,.),  (5.0 

1  1  2  i,j  =  l  1J 

and  let  | V  >  be  the  exact  wave  vector  and  E  the  total 

energy  of  this  system,  i.e. 

H  |  Y>  =  E | Y>  (5.2) 

( 27 ) 

It  can  be  shown  by  variational  method  that  for 

Brueckner  theory  the  exact  wave  function  satisfies  the 
following  equation: 

|  y  >  =  |  <|>  >  +  V  |  'F  >  (5-3) 

where  |  4>  >  is  plane  wave  vector.  In  terms  of  partial  wave 
in  angular  momentum  space,  (5*3)  can  be  written  as 

|  ^(K)  >  =|  ( K)  >  +  — V|  Tt(K)>  (5.4) 

Equation  (5.^1)  can  be  expanded  in  coordinate  space  in  the 
following  way: 


y  0  (K)  >  =  |  4>  (K)  >f  (dt*  dr’  dr"  |K’><K’ 

X/ 


Q 


or  in  wave  function  notation: 


r  co 


y£  (r)  =  <J>£  (r)  + 


TT  JO 


2  ,  Q( K 1 ) 

K  dK  e(K’) 


r ’ > <r  ’  | V  I r^<r" | v|K)> 

(5.5) 

o  ( K  f,r )  4>  0  ( K T  ,  r  ’  )  V  ( r  ’  )  V  ( r  1  )r  ’  ^dr  ’ 


=  <l>£(r)  + 


G£ (r ,r 1  )V(r  ’  )r'H/£(r'  )dr  ’  . 


(5.6) 


.  s  i  i  re  . .  's  .  t  lo  va^ne 


d  -  d 


» 
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Here 


G ’ (r ,r  T )  = 


IT 


<f>  (Kir)<f>  (KT,r')K,2r'  ^  dKT 

o  e  (K*  ) 


(5.7) 


Put  U£(r)  =  rv£(r) 

Then  (5-6)  and  (5-7)  become 


U  (r)  =  r  j  ( Kr )  +  G(r,r')V(r')  U.  (r')dr’ 


(5.8) 


and 


Gn(r>r,)  =  — 


f  rj  (K'r)r ' j  (K'r ' )K'2  Q  (K' }  dK'  (5. 
Jo  e  ( K  ’  ) 


9) 


(K '  ) 

For  potentials  with  hard  core  (5.8)  can  be 
written  as*. 

U£(r)  =  r  j  £  ( Kr )  +  A  G£(r3c)  +j"  G£  (r  3r  ’  )V(r  T  )U£  (r  ’  )dr  » 


with  the  boundary  condition  U£(c)  =  0  where  c  is  the  radius 
of  the  hard  core.  This  implies  that 

cj  £  ( Kc  )  +  AG[c3c)  +1  G£ (r 3r ' )V(r ' )  U£  (r')dr'  =  0 

Jc 


By  solving  for  A 
into  (5.10)  we  get 


from  (5-11)  and  substituting 


Here 


Vr) 


fz(r) 


f£(r)  +  |  K£(r ,r ’ )V(r ’ )U£(r ’ )dr * 


r  j  £  ( Kr )  -  cj  £(Kc  ) 


G  £ ( c ,r ) 

G£  (c  3c  ) 


K&’( r  5  r  '  )  =  G£(r3r  '  )  - 


G  £ (r  3  c )  G£(c3r ’ ) 


G£(c,c) 


.10) 

.11) 

(5.12) 

(5.13) 


and 


<  5  - 1^  ) 


•'  3  >  r.r^liw 


.  v  usp  n  f  rio 


I'-.  -  1  I  •  ,  •  -  )V  ,  ^  (0,3)EM  +  (o«)  to 


•->;  ( 01 . 0  t  (  )ni 
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The  reaction  smplitude  g^  is  defined  by 


s£  =  <f|  V  |Ut  > 

(5.15) 

To  find  the  exact  solution  |U  >  we 

take  a  trial  function 

U£T>and  put 

A[U,t]  -<  St  |v  |uiT  > 

(5.16) 

b[At]  ■  <0tTlv  -  VK  v 

IAt> 

D — 

i — 1 

m 

where  K  V  U£T(r)  =  K£(r,r’)  V  (r’ 

v  0 

)  U£T  (r 1 )  dr' 

(5-18) 

Consider  the  functional 

J[U4T(r)]  = 

CM 

< 

(5-19) 

B 

which  depends  on  the  function  U(r).  As  an  obvious  con¬ 
sequence  of  (5.12  )  and  (5.15  )  we  have 

J[U£(r)]  =  g£  (5.20) 

Here  J[U  (r)]  depends  neither  on  the  normalization  of  the 
function  U  (it  does  not  change  when  U  is  multiplied  by  an 
arbitary  constant),  nor  on  the  value  taken  by  U  in  regions 
V(r)  =  0.  Therefore  J[U  ]  takes  the  same  value  for  any 
function  U  satisfying  the  following  less  restrictive 

X/ 

condition  than  (5*12) 

U£(r)  =  D j  ^ ( Kr )  +  K  V  U£(r)  if  V(r)  t  0  (5.21) 


here  D  is  an  arbitary  constant. 

Now  let  us  calculate  the  variation  of  6  J  as  a 


function  of  6U£T.  From  definition  (5-19) 


6  J  = 


2A 


B 


6  A  - 


A' 


6  B 


(5.22) 


B‘ 


:  iu'  -■  '  •  nr  .  a  r.qs .  ■  -  3.  ri  : 


' -i ;  U'  o  ••  3 


•  jL  .rl  3  '  :  O': 


ta  s  (i^u 


Since  V  is  real  and  K  (r,r')  is  symmetrical  in  r  and  r* 


we 


have , 

6  A  =  <  fj  V  |  6  U£  T  >  (5.23) 

6  B  =  2  <  U£T  |V  -  VK  y  |6  U£T  >  (5.24) 


By  substituting  (5.23)  and  (5.24)  into  (5.22)  we  get 

2A 

SJ  =  —  {  <  ft|V|-S  U£T  ><  UtT  |V  -  VK  V  |UtT>  -<  fj  V  |U1T 

<U4T  |V  -  VK  V  | 6  U(T  >} 

2A 

=  —  {  <  V  lV  -  VK  Vl  °tT  ?  <f*  fJVl  U*T  > 

<  UtTl(V  -  VK  V) | }  S  UlT  > 


2A 


B 


—  (  B  <f  JV  -  A  <U£T  |  (V  -  VK  V)} 


6U£T  > 


2A 

B2 


6  U£T  > 


(5.25) 


where 

<P,,  I  =  B  <fJV  -  A  <U£T|  (V  -  VK  V). 


(5.26) 


In  order  to  have  6J  =  0  for  any  |  6U^T  >,  it  is  necessary 
and  sufficient  that  <  F  | should  vanish.  For  this  it  is 
necessary  and  sufficient  that  <f£|  and  <U£T|  (1  -  K  V) 

be  proportional  at  any  point  where  V(r)  is  not  zero,  that 
is,  be  one  of  the  functions  obeying  equation  (5-21). 

Thus,  the  stationary  value  of  J  is  equal  to  the  reaction 


" )  i 


(  .  )  .‘rtf  (•:.?)  _  U8 


'  WUi  '  •  >  -  ...jU  >■  frfVlj  I  >  -  =  U 
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amplitude  that  we  wish  to  calculate: 

g£  =  extremum  of  J 

If  the  potentials  are  also  dependent  on  spin  and  isospin 
then  the  reaction  matrix  is  related  to  the  reaction 
amplitude  as  follows: 

<K  |  G  |  K  >  =  —5—  H  (2J+1)(J  +  1)  gf’f  (5.28) 

J,T 


and  the  potenteial  energy  per  particle  contributing  to 
the  binding  energy  is 

(5.29) 


E  Kf3 

=  — - -  ZL  (2J  +  1)(2T+1) 

A  6  7t  J ,  T 


K  3  _ 

— —  21  (  2J  +  1 )  (  2T+1 ) 

6  7T  J  J  T  £,d 


The  function  J  can  be  written  in  the  integral  form  in  the 
following  way: 


[|  u*T(r)  v?;j  (r)  ffr)  dr  ]2 

J[U£T]  =  —2 - 

[[  unT(r)  v*;J (r)  u^(r)  dr]  -t  [  u£T(r)v£;J(r) 

Jc  5  Jc  Jc  5 

K£(r3r’  )V^J  (r  ’  )U£T(r  '  )dr ’dr 

(5.30) 


where  f  (r),  G  (r,r')  and  K£(r,r’)  are  given  in  (5-9),  (5.13) 


ST 


and  (5.14). 


:  i  i  2.  j  iR  ic  ob\  i  on  nsffJ 


£  . 


(  '  '   J  -s' 


,  i -;JU  -  [ti-  («»v  <■>■  '  ;>  <*>T, 
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We  choose  the  trial  function  as  the  difference 
between  the  plane  wave  and  a  term  which  corresponds  to 
the  distortion, 


u£T(r)  =  rJ£(Kr) 


cj  z (Kc ) 


h3)(ixtc) 


(5.3D 


which  satisfies  the  boundary  conditions 

U£T(c)  =  0  (5.32) 

U£T(oo)  =  r-MKr)  (5.33) 

In  (5*31)  is  the  parameter  which  dependent  on 
angular  momentum,  spin,  and  isospin. 

In  Pig.  4  the  value  of  J  as  a  function  of  the 
parameter  for  1Sq  is  plotted.  The  stationary  value  of 

J  is  found  to  be  at  the  point  x  =  .4.  Then  from  (5.29) 
we  get  the  contribution  potential  energy  -15.45  MeV  for 

state.  The  value  calculated  from  the  separation  method 
stated  in  previous  chapters  is  -  14  MeV. 


•  '  ;t>I  ■  ;  (  •  '  -,1' 


64 


C-H 


I 


I 


• 

• 

• 

• 

• 

• 

• 

Ch 

Ch 

CJV 

ov 

Ch 

ov 

ov 

VO 

vo 

00 

00 

00 

00 

00 

ro 

o 

00 

Ch 

4^ 

ro 

O 

65. 


Chapter  6. 

CONCLUSION 

Our  results  seem  better  than  others,  (except 
for  the  original  work  of  Brueckner  and  Gammel),  but  the 
agreement  between  theory  and  accepted  value  for  the  binding 
energy  of  nuclear  matter  does  not  imply  the  correctness 
of  the  theory.  For  nuclear  matter,  we  still  lack  enough 
information  to  develop  a  complete  theory.  The  following 
are  some  of  the  fashionable  questions. 

1)  The  uncertainty  about  the  actual  binding  energy  of 

nuclear  matter  because  of  the  extrapolation  of  the  Bethe- 
Weizsdcker  formula  (Green  in  the  International  conference 
on  Nuclear  Phys ics  )  ^  . 

2)  The  region  of  the  validity  of  potential  model  of 
nuclear  forces  -  We  know  that  at  low  energies  the  inter¬ 
action  may  be  represented  by  a  potential,  however,  at 
higher  energies  this  potential  model  may  not  be  valid. 

3)  The  role  of  the  three  body  forces  in  nuclear  matter  - 

We  ignore  these  forces  in  the  theory,  but  lots  of  people 

have  calculated  a  specific  3  body  force  and  the  results 

always  turn  out  to  be  large.  (G.E.  Brown,  International 

(SI) 

Conference  on  Nuclear  Physics  v  Lecture  on  Theory  of 

(16)  j 


Nuclear  Matter 


' 


4)  The  convergence  (or  divergence)  of  perturbation  theory  - 
Here  we  have  assumed  a  rapid  convergence,  but  as  Baker^^ 
has  shown, in  general,  perturbation  theory  diverges.  In  many 
cases,  the  correct  result  may  be  obtained  in  spite  of  the 
divergence  of  the  series.  But  this  cut-off  process  con¬ 
tradicts  the  internal  consistency  of  the  theory  as  in 
quantum  electrodynamics. 

5)  Higher  order  clusters  -  Bethe^^  has  pointed  out  that 
the  three-body  correlation  is  important  and  fourth  order 
is  small  enough  to  be  negligible,  but  the  convergence  (or 
divergence)  for  higher  order  cluster  is  still  not  clear. 

Although  the  theory  of  nuclear  matter  is  still  far 
from  being  comprehensive,  we  agree  with  Bethe  that  it  is 
at  last  at  a  point  where  "it  can  be  used  with  some  con¬ 
fidence  to  calculate  quantities  which  have  not  been  accurately 
observed,  like  compressibility,  symmetry  and  surface  energy, 
and  where  it  is  also  applied  to  finite  nuclei." 
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